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We develop a general Shockley-Ramo-type framework to describe spatial patterns of photocurrent 
response in gapless materials. Our approach helps to understand the striking features of the observed 
patterns, such as the directional effect and the global character of photoresponse. We illustrate this 
approach by examining specific examples, and show that the photoresponse patterns can serve as a 
powerful tool to extract information about symmetry breaking, inhomogeneity, chirality, and other 
local characteristics of the system. Short response times, originating from the nonlocality, make 
photocurrent uniquely sensitive to the charge and spin dynamics of the system. 



The optoelectronic response of many gapped materi- 
als features spatial nonlocality arising due to slow re- 
combination of photoexcited carriers [lj. Recently, scan- 
ning photoresponse techniques were used to probe new 
gapless materials, such as graphene and topological in- 
sulators, where carrier recombination is fast compared 
to the characteristic timescales for transport across the 
system (sec Refs.jH-S] and references therein). Surpris- 
ingly, these studies revealed complex spatial patterns of 
photocurrent response. Rather than being localized near 
current-collecting contacts, the photocurrent hot spots 
are distributed over the entire system area, typically sep- 
arated by many microns from the contacts [3|-|7|. These 
large length scales may seem hard to reconcile with the 
short picosecond-scale recombination times over which 
the photoexcited carriers lose their energy and become 
part of the thermal distribution, traversing distances 
much less than the system size. 

Here, we introduce a framework that naturally explains 
how the non-locality can arise in the absence of slow 
recombination, and incorporates diverse local phenom- 
ena predicted theoretically Using this framework, 
we can understand striking features of the observed pat- 
terns 0-0, such as the global character of the response 
and the directional effect (Fig[T}. Namely, the photocur- 
rent hot spots are highly sensitive to the orientation of 
inhomogeneitics and interfaces, at which the hot spots 
are pinned, but totally independent of the distance from 
the contacts [FigOJd) shows data from Ref. 0] where this 
effect was first reported]. Further, the resulting pat- 
terns are sensitive to the symmetries which govern pho- 
toresponse. This is illustrated for several mechanisms, 
charge-dominated (FigJlJ and spin-dominated (Fig(3j) . In 
the latter case, we show that the spatial pattern structure 
is directly sensitive to the details of spin dynamics. 

In gapless materials, photoresponse is mediated by am- 
bient carriers outside the photoexcitation spot. Short 
recombination times lead to a rapid decay of the pri- 
mary photoexcited carriers, preventing them from reach- 
ing contacts and directly contributing to photocurrent. 
Instead, the main contribution is indirect: a local pho- 
tocurrent sets up an e.m.f. that drives ambient carriers 




FIG. 1: Scanning photocurrent images for different mecha- 
nisms of photoresponse. The photocurrent, drawn from con- 
tacts 1 and 2, is modeled by Eqs.([T]),([5}. (a) Photocurrent 
pattern in a chiral material, where an = ±1 marks regions of 
different chirality. Local photocurrent direction is governed 
by edge states (white arrows), (b) Photocurrent pattern in 
a non-chiral system with a step-like density inhomogeneity 
(see text), (c) Photocurrent pattern due to floating contacts 
that do not draw current (yellow semi-circles labeled 3-8). In 
all 3 cases, the system-wide photoresponse has strong direc- 
tional dependence on the local photocurrent orientation and 
profile, but essentially no dependence on its location within 
the system, (d) Scanning photocurrent image of a 12 /im-long 
graphene device with six floating contacts 3-8. Note that the 
sign of photoresponse near floating contacts is correlated with 
the direction to the contacts 1 and 2 through which current 
is drawn (data taken from Fig.2(a) of Ref. [J]). 

far from the excitation region, and into the contacts, lead- 
ing to response which does not diminish with distance. 

The "global" character of the response can be linked 
to charge continuity. The condition V • j = can be inter- 
preted as incompressibility of current flow, with stream 
lines that do not terminate anywhere within the system. 
In addition, because the current is caused by a chemical 
potential gradient, the stream lines cannot close on itself. 
This results in a response at contacts far away from the 
spatial structures probed by the photoexcitation. As we 
shall see, the system-wide response to local photoexcita- 
tion can be described as 

I = A /j loc (r)-W'(r)d 2 r, (1) 
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where ji oc ( r ) is local photocurrent in the photoexcita- 
tion region, ip is a weighting field obtained by solving a 
suitable Laplace problem, and A is a prefactor which de- 
pends on device configuration (see Eq.@). Spatial pat- 
terns predicted using Eq.([T]) exhibit photocurrent- active 
structures with contrast which is essentially independent 
on their position within the system (see FigUJ. Such 
"global" photoresponse is known for one-dimensional sys- 
tems, where Eq.([IJ reduces to adding up the total poten- 
tial drop across the device However, the generalized 
framework presented here yields photocurrent that can 
exhibit complex structures which are not anticipated in 
a one-dimensional approach. 

There is an interesting relation between our result and 
the seminal Shocklcy-Ramo (SR) theorem developed in 
the context of vacuum-tube electronics ljl li| . Shockley 
and Ramo were concerned with the currents induced in 
the electrodes by charges moving in the free space inside 
the tube. The SR relation, which is intrinsically nonlocal 
due to the long range character of electric field in vacuum, 
can be expressed in a closed form as 

I jk oc ev(t) ■ E w | r=R ( t ), E„(r) = Vto jfc (r), (2) 

where Ijk is the current induced by a moving charge e 
in the electrodes j and k, where v(t) and R(i) are the 
charge velocity and position. The function Wjk(j) sat- 
isfies Laplace's equation with suitable boundary condi- 
tions on the electrodes. The SR theorem, which can be 
extended to charges moving in insulators [1 71] . is key in 
ultra-fast charge sensing, such as particle detection in 
high energy physics flil plasma diagnostics po}. etc. 

We emphasize that the origin of nonlocality in our pho- 
toresponse problem is quite different from that in the SR 
problem, since the ambient carriers screen the long range 
electric field created by photoexcited carriers. Instead, 
the nonlocality arises due to long range currents. Fur- 
thermore, the SR theorem is typically applied to high- 
speed charge detection, whereas we are concerned with 
the steady-state photocurrent. Yet, despite these dif- 
ferences, our approach yields a relation [Eq.([T])] which 
exhibits formal similarity with the SR theorem. 

The starting point of our analysis is the continuity 
equation, V ■ (jd + ji oc ) = 0, where is the diffusion 
current due to ambient carriers in the material. These 
two contributions to current have very different spatial 
dependence: the photocurrent ji oc is present in the exci- 
tation region, whereas the diffusion current jd is nonzero 
throughout the entire material, describing the response 
to an c.m.f. due to ji oc driving current far outside the 
excitation region. Below we focus on the simplest sit- 
uation when transport can be described by a position- 
dependent conductivity tensor a(r) (other effects, such 
as the hot carrier transport [2l|, [22J, can be analyzed 
similarly). The diffusion current satisfies the usual re- 
lation ]d = — c(r)V0, where (f> is the electrochemical 
potential. The boundary conditions in this transport 



problem are zero current through the sample boundary, 
n " (jd +jioc) = 0, and constant potential at the contacts, 
n x = (we assume a two-dimensional geometry). 

To handle the non-local response, we introduce an aux- 
iliary weighting field ip(r) satisfying V • j^ (r) = in the 
bulk of the material, where jw = — a T Vip, and appropri- 
ate boundary conditions at the boundary and contacts, 
n . jW(r) = and n x W( r ) = 0, respectively. Mul- 
tiplying the continuity equation for the physical current 
jd + jioc by ip{ r ), integrating over the sample area, and 
using Gauss' theorem, we obtain 



W(r) -jioc(r)d 2 r 



k 



(3) 



where k labels contacts. The quantities on the right hand 
side are the net currents flowing in each of the contacts, 
Ik = J c n • jkd£, and potentials on these contacts. We 
emphasize that Eq.(|3]) holds on very general grounds re- 
gardless of whether a particular contact is drawing cur- 
rent (Ik ^ 0) or is floating (Ik = 0). The expression on 
the left hand side depends on the microscopic photocur- 
rent distribution ji oc (r) inside the material, whereas the 
expression on the right hand side is a function of currents 
and potentials at the contacts, thereby providing a gen- 
eral relation between position-dependent photoexcitation 
and the measured photocurrent. 

It is convenient to choose ip(r) such that 1^ = for 
all floating contacts. Then the contribution to Eq.Q due 
to floating contacts drops out entirely, yielding a relation 
which only includes the contacts that actually draw cur- 
rent. It is also straightforward to account for the effect 
of an external circuit. We consider the current drawn 
through a pair of contacts 1 and 2 (see Fig[TJ and write 

4p) = ~ <h)/R, h(2) = ±(0i - <h)/R**t, with R 
and i? cx t the resistance of the sample and of the exter- 
nal circuit, respectively. Setting tpi — ip2 = lj we obtain 
Eq.(p} with the prefactor 



A = R/(R + R ext ). 



(4) 



Despite its apparent simplicity, Eq.(UJ accounts for the 
key effects that impact photoresponse, such as system 
geometry, structure, inhomogeneity, chirality, etc. 

The general features of Eq. flJJ can be illustrated for a 
spatially uniform system of a rectangular shape. In this 
case, the weighting field ^>(r) is a linear function, \7ip = 
y/L, with L the system length (see FigJT]). Constant 
Wi/j leads to a "global" response which is invariant upon 
spatial translation of ji oc ( r )- At the same time, the sign 
and the magnitude of the response depend on the angle 
between Vip(r) and ji oc (r) (the directional effect). 

The value ji oc (r) depends on system properties in the 
photoexcitation region. By symmetry, no photocurrent 
can occur in a spatially uniform system (assuming un- 
polarizcd light). In the presence of a density gradient 



Vn(r), the local photocurrent can be described by 
ji oc (r) = [aVn(r) Vn(r)] J(r), 



(5) 



where a and /3 are material constants, and J(r) is the 
absorbed optical power. In general, a is finite in all ma- 
terials, whereas /3 is only non-zero in chiral systems where 
edge-state transport allows ji oc to be directed along the 
contours of n(r). This is the case in chiral materials such 
as topological insulators due to coupling between orbital 
motion and spin [H, [ll , 12 1 , or in non-chiral materials in 
the presence of a magnetic field 

The effects of spatial inhomogeneity arc illustrated in 
FigH] for the chiral response (a) and the nonchiral re- 
sponse (b). For illustrative purposes, we use a step-like 
density profile, with n taking one value in the middle re- 
gion and another value in the top and bottom regions, 
identical for (a) and (b). For (a) we use j loc with a = 
and finite /3, for (b) it is the other way around. In both 
cases, the photocurrent is zero in the regions of constant 
n, and nonzero near the steps. The differences in the sign 
and magnitude of the response reflect the fundamental 
difference in physics in the cases (a) and (b). 

Model (a) describes photoresponse in a chiral systems, 
which arises at the interfaces between domains of op- 
posite chirality, labeled by <tjj = ±1. Physically, it 
may represent a quantum Hall system near a plateau 
transition[23j], or a system in which nonzero chirality re- 
sults from spontaneous ordering The different signs 
of chirality can be associated with the clockwise and 
counter-clockwise edge states, labeled by white arrows. 
Notably, the sign and magnitude of photocurrent depend 
on the direction of current flow in the edge states. The 
photocurrent is also nonzero at system boundaries, indi- 
cating the presence of current carrying edge states. This 
can be used to identify the edge states and domains with 
different chirality in experiment. 

FigfT] (b) shows the non-chiral photocurrent response 
for the same density profile as in FigfT] (a). Physically, 
(b) may describe systems such as graphene with spatial 
inhomogeneity giving rise to p-n boundaries separating 
regions with electron- like and hole-like polarity In 
this case, ji oc is normal to the contours of n(r), making 
the sign and magnitude of the response dependent on the 
orientation of the interfaces viz. y - jioc- Also, since ji oc 
is normal to boundaries whereas VV> is tangential, the 
photocurrent vanishes at the system edge. 

A very different behavior is found near contacts, since 
VV> is normal to the contact surface, see FigQ] (c). In 
this case, a nonzero response arises both near the con- 
tacts through which current is drawn and near floating 
contacts (see also FigfJ]). Notably, the response depends 
on the floating contact orientation but not on its position 
within the system. This is in agreement with experimen- 
tal observations of Ref. [1] , which are reproduced in Fig. 
[TJ (d) . All three photocurrent patterns in FigfT] despite 
their different physical origin, share two common trends: 
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FIG. 2: Directional effect in photoresponse accounting fully 
for the distortions of the weighting field, (a) Photocurrent 
pattern due to three circular regions, modeled in the same 
way as in FigfTJ (b). The conductivity inside each region is 
taken to be 10 times larger than the background conductivity. 
(b,c) Photoresponse and the field lines for Vt/> near floating 
contacts of two different shapes, a rectangle and a semicircle, 
obtained using the conformal mapping approach, Eq.©. 

strong directional sensitivity and global character (posi- 
tional independence). This behavior makes the photocur- 
rent patterns particularly useful in identifying symmetry 
breaking and inhomogeneity in gapless materials. 

The patterns in FigfTJ (a), (b) were obtained using a spa- 
tially uniform weighting field approximation, \7ijj « y/L. 
To analyze a more realistic inhomogeneous situation, we 
used a numerical procedure to obtain the exact profile 
-0(r). Figf2] (a) shows photocurrent patterns from three 
circular regions with a mismatch between the inner and 
outer conductivity, which causes significant distortions of 
the Vip field lines. Yet these distortions do not impact 
the overall trends discussed above, the global character 
of the response and the directional effect. 

In contrast, the weighting field distortions have a very 
dramatic effect near contacts. Even if a contact does not 
draw net current, it short-circuits the current flowing in 
its vicinity, leading to a non- vanishing normal component 
of V'0 near the surface of a contact (see Fig. f2]). For ji oc 
which is normal to the contact, this gives a nonzero, sign- 
changing photoresponse, as in FigfT] (c,d). 

For ideal contacts, the field ip can be found using the 
conformal mapping approach, giving ^(r) = Almw(z). 
Here w is a suitable analytic function of a complex vari- 
able z = x+iy, which satisfies the equipotential condition 
at the contact surface. We illustrate this for a flat contact 
and for a semicircular contact (see Figf5] (b,c)): 



(z) = V( z ~ Vi)( z - D2), w c (z) = z- Zq 



- zo 
(6) 

where the flat contact is positioned at yi < y < y%, x = 0, 
and the semicircular contact is of radius r and is posi- 
tioned at z = zo. Here we assume that the contacts 
are floating and are small compared to the system size. 
At large z, ip asymptotically approaches the linear de- 
pendence ip oc y found above. The photocurrent at the 
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FIG. 3: Spatial dependence of spin-driven photocurrent, 
Eq.©,©, for spin diffusion length values £ = 5L, L/10, L/25 
(curves 1,2,3). Spin-diffusion length £ can be directly deter- 
mined from the measured photocurrent profile. 

contact is proportional to n • V-0. For the flat contact, 



A- 



y- \{vi + m) 
V(y-yi)(y2-y)' 



yi < y < V2- (7) 



Since this quantity is an odd function of y— h (y± +2/2) , the 
net current drawn the contact vanishes, as appropriate 
for a floating contact. Similar sign-changing behavior is 
found for the semicircular contact, see Fig [2] (c). The 
sign-changing pattern is oriented in such a way that the 
parts with high photoresponse are facing the contacts 
through which the photocurrent is drawn, in agreement 
with the observations reported in in Ref. j4j , see Figfjjd) • 
The approach described above can be applied to a 
variety of photocurrent mechanisms. As an illustra- 
tion, we discuss the spin-driven photocurrent which arises 
in materials without inversion symmetry, in particular 
the topological insulators @, [H, (2(| and zinc-blende 
systems [27]]. In such systems, spin orientation by circu- 
lar polarized light results in spin imbalance S(r) oc J(r), 
generating local photocurrent via 



j loc (r) = QS(r). 



(8) 



Here Q is a second-rank tensor describing the spin- 
galvanic effect. Since the system-wide photoresponse, 
Eq. ([!]). is sensitive to spatial distribution of spin, it can 
be used to probe spin dynamics and relaxation. 

As an illustration, we consider the rectangular geome- 
try of FigQ] and assume weak spin flip rate. In the diffu- 
sive limit, spin transport is governed by -DV 2 S + S/t = 
>^s [Hi j where Js describes spin orientation by absorbed 
light. Evaluating the photocurrent, Eq.([T|), we find 

= (cosh A _ cosh / ( C0S h L _ A (9) 

where £ = \J Dr. Notably, for £ <C L, the current is 
finite even when spins decay before reaching contacts. In 
contrast to the charge-dominated response, photocurrent 
vanishes near contacts which act as a spin sink. The 



measured profile I(y) can be used to directly determine 
the spin diffusion length (see Figj3]). 

Importantly, ambient carriers lead to a very fast pho- 
toresponse. A simple estimate of timescales can be given 
by reinstating the time dependent term in the continuity 
equation. For a spatially uniform system, the dynam- 
ics can be expressed in terms of the Fourier harmonics of 
charge density, d t 8n^{t) = — 2^er|k|<5nk(i), where a is the 
sheet conductivity per square area and k is the dielectric 
constant. For a crude estimate, taking parameter values 
|k| ps tt/L, L = 10 fim, k = 5, 1/cr = 1 kfi, we obtain a 
sub-picosccond response time 



t = kL/(2tt 2 <j) ps 0.3 ps, 



(10) 



which is considerably shorter than typical cooling and re- 
combination times. Fast response makes the system-wide 
photocurrent a useful probe for the dynamical processes 
in the excitation region. It also makes gapless materials 
viable for applications in high-speed optoelectronics. 

In summary, our approach explains several puzzling 
aspects of photocurrent response in gapless materials, 
in particular the striking non-locality and the direc- 
tional effect observed in Ref.Q. By analyzing different 
mechanisms of photoresponse, we demonstrate that it 
is uniquely capable of revealing spatial patterns arising 
due to symmetry breaking, chirality, or inhomogeneities. 
Sub-picosecond response times make photoresponse a 
useful probe of system dynamics. 
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